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ABSTRACT 

Context. The most massive stars are thought to be hydrogen-rich Wolf-Rayet stars of late spectral subtype (in the following WNh 
stars). The emission-line spectra of these stars are indicative of strong mass loss. In previous theoretical studies this enhanced mass 
loss has been attributed to their proximity to the Eddington limit. 

Aims. We investigate observed trends in the mass-loss properties of such young, very massive stars to examine a potential In- 
dependence, i.e., with respect to the classical Eddington factor T c . Based on different mass estimates, we gain information about 
the evolutionary status of these objects. 

Methods. We derive theoretical mass-luminosity relations for very massive stars, based on a large grid of stellar structure models. 
Using these relations, we estimate Eddington factors (T e ) for a sample of stars, under different assumptions of their evolutionary 
status. We evaluate the resulting mass-loss relations, and compare them with theoretical predictions. 

Results. We find observational evidence that the mass loss in the WR regime is dominated by the Eddington parameter T c , which 
has important consequences for the way we understand Wolf-Rayet stars and their mass loss. In addition, we derive wind masses that 
support the picture that the WNh stars in young stellar clusters are very massive, hydrogen-burning stars. 

Conclusions. Our findings suggest that the proximity to the Eddington limit is the physical reason for the onset of Wolf-Rayet type 
mass loss. This means that, e.g. in stellar evolution models, the Wolf-Rayet stage should be identified by large Eddington parameters, 
instead of a helium-enriched surface composition. The latter is most likely only a consequence of strong mass loss, in combination 
with internal mixing. For very massive stars, the enhanced Independent mass loss is responsible for the formation of late WNh sub- 
types with high hydrogen surface abundances, partly close to solar. Because mass loss dominates the evolution of very massive stars, 
we expect a strong impact of this effect on their end products, in particular on the potential formation of black holes, and Gamma-Ray 
Bursts, as well as the observed upper mass limit of stars. 

Key words. Stars: Wolf-Rayet - Stars: early-type - Stars: atmospheres - Stars: mass-loss - Stars: winds, outflows 

1. Introduction pends on the Eddington factor F e (Eq.Q]). Such relations have 

been predicted for very massive s tars close to the Eddington 



Wolf-Rayet (WR) type mass loss fundamentally affects the limit jGrafener & Hamannl l2008t IVink et all IZOlTh . and for 



evolution, the final fate, and the chemical yields of massive LBVs dVink & de Koterf l2002h . Notably, the proximity to the 



stars. The amount of mass loss in the WR phase predomi- Eddington limit provides a natural explanation for the occur- 

nantly decides whether a star ends its life as neutron star, or rence of the WR phenomenon. Large Eddington factors can be 

black hole dHegeret al.| |2003|). In particular, WR-type mass reac hed, on the one hand, by very massive stars on the main se- 

03 ; loss at low metallicities (Z) is expected to be of paramount quen ce because of their extremely high luminosities, and on the 

importance for the chemical enrichment of the early universe other handi by less mass i ve evolved (He-burning) stars due to the 

dMeynetetalJ |2006t |Chiappini et alj |2006|), and the formation enhanced mean molecular weight in their cores. The occurrence 

of long-duration gamma ray bu rsts (long GRBs, |Yoon & Langerj of W R-type mass loss for young, luminous, hydrogen-rich stars 



|2005jlWoosley & HegeiEQQj). (typically late WNh subtypes), and evolved, hydrogen-free WR 



The nature of WR-type mass loss is however still poorly stars (nitrogen-rich WN, and carbon-rich WC subtypes) can thus 

understood. Stellar evol ution models mostly rely on empirical be explained in the same way. 
mass- loss relations (e.g. iNugis & Lamersl |2000T : lHamann et al.l 

l2006h. with the WR phase identified on the basis of observed The main goal of this paper is to confront the theoretically 

WR surface abundances in our galaxy. Clearly, for the model- predicted concept of T-dependent mass loss with the observed 

ing of stellar populations that cannot be observed locally, a more mass-loss properties of very massive stars. We use the results of 

physical approach would be desirable. a stu dy of the most mas sive stars in th e Arch es cluster near the 

In the present work we elaborate on such an approach, Galactic centre (GO. bv lMartins et al] d2008l) . The Arches clus- 

namely a mass-loss relation for WR stars that chiefly de- ter is a young star forming region rich in very massive stars, in- 
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eluding many O and Of+ supergiants, and luminous WNh stars. 
It forms an ideal testbed for our present study. 

Our paper is organized in the following way. In Sect.[2]we 
briefly describe the properties of luminous WNh stars. Sect. [3] 
recaps our theoretical knowledge of the mass-loss properties of 
these objects. To study the dependence of mass loss on T e re- 
quires mass estimates for the sample stars. In Sect.|4]we provide 
mass versus luminosity relations for chemically-homogeneous 
stars, that may be used to derive the Eddington factor for the 
Arches stars from their observed luminosities, and surface abun- 
dances. In Sect. [5] we study the dependence of the mass-loss 
rates on stellar properties, including T e , and compare these 
dependencies with theoretical predictions. The main intention 
in this section is to perform a qualitative comparison. At the 
present time a quantitative comparison may be affected by sys- 
tematic uncertainties. Our findings are discussed in Sect.|6l and 
conclusions are drawn in Sect.|7] 

2. Physical properties of the most massive stars 

The most massive stars with direct mass estimates are found 
in binary systems that contain luminous, hydrogen-rich WNh 
stars. The highest masses lie in the range of 70-120 Mr 



stars, ine nighest masses lie in the range ot /u-lzUM a 
dRauw et alJll996HSchweickhardt et al.||1999tlRauw et al.ll2004t 



iBonanos et all 12004 ISchnurr et al.l 120081 l2009h . Furthermore 
spectral analyses of (putatively) single WNh stars imply very 
high masses. Typical luminosities lie in the range of 1O 6 L 
or higher, implying that most WNh stars are very mas- 
sive st ars in the phase of cor e H-burning dCrowfher et aT 



19951; ICrowther & Smifhl 119971: ICrowther & Dessartl 



de Koter etaD 11997b iHamann et all 120061: iMartins et ail 
20091 ICrowther et alJl2010l) . 



1998 



2008 



Luminous WNh stars are preferentially found in the cen- 
ters of massive, young stellar clusters with ages of only a 
few Myr. This, combined with the fact that their surfaces still 
show ample amounts of hydrogen, suggests that they are still 
in their phase of core H-burning. Well-known examples are 
the Arches cluster close to the Galactic centre, the young 
galactic cluster NGC3603, and R136, the central cluster of 
the star -forming region 30 D or in the large Magellanic cloud 
(LMC). ICrowther et al] d20 1 0b recently determined a luminosity 
of 1O 694 L for the brightest star in R 136, which would corre- 
spond to a single star with a present-day mass of 265 M Q . 

The fact that the most massive stars appear as Wolf-Rayet 
spectral types can be explained as a result of their prox i mity t o 
the Eddington limit. According to IGrafener & Hamannl d2008h . 
the increased density scale height close to the Eddington limit 
leads to the formation of strong winds with large optical depths. 
Because ionizing photons are efficiently absorbed within these 
winds, recombination sets in. This shows up in the form of strong 
WR emission lines that originate from the subsequent recombi- 
nation cascades. 

In the following we characterize the proximity to the 
Eddington limit by the 'classical' Eddington factor 

T e =xM(4*cGM), (1) 

with the electron scattering opacity x$H Because H, and He are 
completely ionized in the inner regions of hot star atmospheres, 
T e is nearly constant with radius, and depends only on the stellar 
parameters M, L, and, via the mass absorption coefficient Xe, on 
the hydrogen mass fraction at the stellar surface (Eq.[8]l. 



There exists a significant additional contribution to the total 
mean opacity x( r )> due to metal lines (chiefly Fe) and continua. 
For the physical, radius dependent Eddington factor T(r), includ- 
ing all opacities, we thus have 



T{ r )=x(r)LI{4ncGM)>T e . 



(2) 



The size of the shift between T and T e depends on the detailed 
metal abundances, and on the ionization structure of the atmo- 
sphere, and thus on r», and Z. E.g., for solar metallicity models, 
IGrafener & Hamannl (120081) find that V approaches unity in deep 
atmospheric layers, already for F e « 0.5. The onset of WR-type 
mass loss thus occurs for T e < 1, i.e., for Eddington factors that 
are smaller, but still of the order of one. 



3. Mass-loss predictions for very massive stars 

Mass-loss predictions for very massi ve sta rs have been per - 
formed by IGrafener & Hamannl d2008l) . and IVink et all d201 ll) . 
Despite the rather different modeling approaches, these works 
agree on the dominant role of the Eddington factor T e for the 
mass-loss properties of very ma ssive stars. 

The mass-loss relation by Gra fener & Hamannl d2008l) is 
based on advanced stellar atmosphere models that incorporate 
non-LTE line blanketing, wind clumping, a nd an exact numerical 
soluti on of the hydrodynamic equations IGrafener & Hamannl 
120051) . The models include complex model atoms of H, He, C, N, 
O, Si, and the Fe-group, which should be sufficient to describe 
the largest part of the radiative wind acceleration. We note how- 
ever that the lack of intermediate elements (Ne-Ca) could poten- 
tially lead to an under-estimation of th e mass-loss rates. 

Apart from the dominant role of r„. IGrafener & Hamannl find 
a strong dependence on T*, and a strong Z-dependence. Their 
mass-loss prescription has the form 

log(M) = -3.763 

+ p ■ log (T e - r () ) - 3.5 ■ (log(7VK) - 4.65) (3) 
+ 0.42 ■ (log(L/L ) - 6.3) - 0.45 ■ (X H - 0.4) , 

with M in M yr , and the Z-dependent parameters /?, and To 
given by 

/3(Z) = 1.727 + 0.250 • log(Z/Z ), (4) 
F (Z) = 0.326 - 0.301 • log(Z/Z ) - 0.045 ■ log(Z/Z ) 2 . (5) 

Note that the stellar temperat ure T+ denotes the eff ective core 
te mperature as define d, e.g., in IGrafener et al.l (l2002h . 

Vink et al.l d201 ll) compute mass-loss rates for very massive 
stars using a Monte Carlo approach with a parameterize d solu- 
tion of the hydrodynamic equations dMiiller & Vink|2008l) . They 
confirm the dominant role of F e , but find a we ak tempera ture de- 
pendence in the range of 30-50 kK. Notably. IVink et ail resolve 
the transition between classical OB star mass loss, with a rela- 
tively weak dependence on F e (Eq.|6j, and WR-typ e mass loss 
with a much steeper dependence (Eq.|7j. However, IVink et aT] 
note that the precise value of T e where this transition occurs in 
their models, might be too high. The likely reason is that the shift 
between the Eddington factor T(r), and F e (cf. Eq.|2|i is under- 
estimated in their models (see also the discussion in Sect. l6.2b . 

F or a solar c omposition, and an effective temperature of 
50 kK lVink et al.1 give relations of the form 



1 Note that the opacity ^ is a mass absorption coefficient, i.e., in the 
CGS system it is measured in cm 2 /g. 



log(M) oc 1 .52 log(F e ) + 0.68 log(L/L ) T e < 0.7 



(6) 



log(M) oc 3.99 log(F e ) + 0.78 log(L/L Q ) T e > 0.7, (7) 
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with M in M Q yr -1 . 

For the Eddington factor F e both works adopt the value for a 
fully ionized plasma, which is given by 

logCTe) = -4.813 + log(l + X S H ) + log(L/L ) - log(M/M Q ). (8) 

In this form F e only depends on the stellar parameters M, L, and, 
via Xe, on the hydrogen mass fraction X H at the stellar surface 
(cf. EqJD. 

The strong sensitivity to the Eddington factor in both mass- 
loss relation^ offers the potential to provide very precise esti- 
mates of r e for specific objects. If M, L, X H , and Z are known 
from spectral analyses, it is thus possible to obtain very precise 
mass estimates, within the systematic errors of the adopted mass- 
loss relation. 

A major goal of the present work is to calibrate the under- 
lying mass-loss relation by a comparison of such 'wind masses' 
with predicted masses from stellar structure computations. If the 
important dependencies on T e , T*, and Z are backed up by ob- 
servations, the mass-loss relations may not only serve as input 
for stellar evolution computations, but can also provide an im- 
portant diagnostic tool to examine the present masses, and thus 
the evolutionary status of observed stars. 

4. Mass-luminosity relations for very massive stars 

In the present section, we provide theoretical M-L relations 
for very massive stars. These relations can be used to estimate 
masses M, and Eddington factors T e for observed stars with 
known stellar parameters L, T+, and X H . 

A basic problem of such an approach is that the internal 
structure of a single observed star is generally not known, i.e., 
its precise mass cannot be uniquely predicted. In the present sec- 
tion we thus focus on the extremes, namely the lowest and high- 
est masses for a star with given observed parameters. Applying 
these relations to a large sample of stars, we will be able to per- 
form a qualitative investigation of their mass-loss properties, and 
to examine patterns that are related to the stellar core, and sur- 
face abundances. 

The highest possible mass for a star with given luminosity L, 
and hydrogen surface abundance X H is reached by chemically- 
homogeneous stars. Under this assumption, the star is character- 
ized by one, constant hydrogen abundance Xu, which equals the 
surface abundance X^. The estimated stellar mass M^omiL, X R ) 
is strongly dependent on the (observable) surface abundance X H . 
Stars that are not homogeneous have a higher mean molecular 
weight in the core than at the surface. They thus have higher 
L/M ratios, or lower masses, for given L, and X^. 

The lowest possible mass is given by the completely inhomo- 
geneous case, i.e., by a core hydrogen abundan ce X^ = 0. In this 
case th e star is in the core He-burning phas^l Lauterborn et alj 
dl97ll) have investigated how luminosities, and temperatures of 
core He-burning stars vary, depending on the mass ratio between 
the H-rich envelope and the He core. According to this work, 
there exists a generalized main sequence that depends on the 
size of the He core. Stars with very small He cores display hot 
temperatures, and similar luminosities as their H-burning coun- 
terparts. For larger cores, temperatures become very cool, and lu- 
minosities increase with increasing core size. Finally, stars with 

2 The steepness of the r c -dependence in Eq.[3]is mainly due to the 
fact that log (r e - r ) is evaluated for T c - T close to zero. 

3 Note that we only focus on stars in the phase of central H, or 
He-burning because the burning timescales in later phases become ex- 
tremely short. 




1.0 1.5 2.0 2.5 

log w (MIMq) 



Fig. 1. Homogeneous stellar structure models for the mass range 
8-250 M (black symbols). The models are computed for hydro- 
gen mass fractions Xu = 0.7-0 (from bottom to top). Dashed red 
lines indicate the M-L relations according to Eq.|9]& row No. 1 
in Table lA. ll for Xu = 0.7, 0.4, and 0.1, and the solid red line cor- 
responds to pure He models according to Eq.[10]& ro w No. 6 in 
Table l A. II For comparison, model computations from llshii et al.l 
(1 19991) are indicated by red circles. The fitting relations are in- 
ferred for the mass ranges of 1 2-250 M Q (Eq.gji, and 8-250 M G 
(Eq.HOjl. 



large He cores show the same luminosities as pure He-stars, and 
hot temperatures. Because we are interested in hot, massive stars 
with large convective cores in this work, we can thus use the M- 
L relation for pure He-stars for the core He-burning case. Our 
minimum mass is thus given by Mueb(L) = Mh om (L, Xu = 0). 
Note that this mass estimate is completely independent of the 
(observed) surface abundance X R . 

The different dependence on X^ in both cases is of 
paramount importance for the present work. Because of the ex- 
pected relation between mass and mass loss, we will be able to 
identify patterns in the observed properties of very massive stars 
that help to distinguish between samples of well-mixed stars, 
that are quasi-chemically homogeneous, and stars with a pro- 
nounced chemical profile. 

4. 1 . The M-L relation for homogeneous stars 

In the present section we derive analytical expressions of the 
form L(M, Xu), and M(L, Xu), for the masses and luminosities of 
chemically-homogeneous stars. As discussed above, we expect 
the luminosity L to depend on the mass M, and the hydrogen 
mass fraction Xu for core H-burning stars, and only on M for 
core He-burning stars. To investigate this dependence we have 
computed a grid of homogeneous stellar structure models for 
the mass range M = 0.3^-000 M G , and hydrogen mass fractions 
Xu = 0.0-0.7, at solar metallicity. 

The stellar structure models are computed with a simple code 
that integrates the stellar structure equations using a shooting 
me thod. The numerica l meth ods are described in the textbook 
by lHansen & Kawaleil dl994l) . The code is based on an exam- 
ple program that is distributed with the book , but is completely 
re-wr itten, and updated with OPAL opacities (llglesias & Rogers] 
Il996l) . To circumvent numerical problems with the extended 
envelopes of extre mely massive stars (see llshii et ail 1 1 999t 
iPetrovic et al.ll200"6h . we have adopted an outer boundary tem- 
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Fig. 2. M-L relation for chemically-homogeneous stars, depen- 
dence on Xu- The symbols are the same as in Fig.Q] Red lines 
indicate our linear fit to models with constant mass but varying 
hydrogen mass fraction Xu, for a mass range of 12-250M o , ac- 
cording to Eq.[9]& row No. 1 in Table lA. ll 



perature that lies above the temperature of the Fe-opacity peak 
(~ 160 kK). In this way we start our computations just below the 
extended envelopes. Because the masses of such envelopes are 
very small (< 10~ 2 M*) this approach has no effect on the ob- 
tained luminosities. 

The results of our grid computations are displayed in 
Fi gs.[TJ&|2l They are i n good agreemen t with previous models 
by Langei d!989l) . and llshii et all ( 1 19991) that comprise a much 
smaller parameter range. We find that a polynomial relation that 
is quadratic in log(M), and linear in Xu fits the results satisfac- 
torily over the parameter range M = 12-250 M , and Xu = 0.1- 
0.7 (see Figs.[TJ&[2]i. The resulting relation has the form 



log(L/Lo) = [Fi +F 2 X H ] 

+ [F 3 + F 4 X H ] log(M/M ) 
+ [F 5 + F 6 X H ] log(M/M ) 2 



(9) 



with the coefficients F\-F(, from row No. 1 in Table lA. ll For 
Xu < 0.1, the changes in the core temperature become so large 
that the dependence on Xu becomes significantly non-linear. We 
thus restrict our fit to Xu > 0.1, and derive a separate relation 
for pure He stars {Xu - 0), over a mass range of 8-250 M Q . This 
relation has the form 

log(L/L ) = F l+ F 2 log(M/M ) + F 3 log(M/M ) 2 , (10) 

where the coefficients F1-F3 are again given in row No. 6 in 
Table lA. ll The maximum fitting error for log(L/L ) over the 
given parameter ranges amounts to 0.02. Inverting these rela- 
tions we obtain the masses M^omiL, Xu), and Mueb(L) in the 
form 



F] +FoX H + F 3 Jf 
log(M hom /M ) = yJ 

1 + tg Xu 

with 

/ = F A + F 5 X n + F 6 X 2 + (F 7 + F & Xu) log(L/L ), 
and 

log(M He b/M ) = F l+ F 2 V^3 + mog(L/L ). 



(11) 



(12) 



(13) 



All coefficients are given in Table lA. ll together with coefficients 
for additional relations that cover higher, and lower mass ranges. 
The latter are discussed in Appendix [Al and may be useful for 
future studies. 



5. The most massive stars in the Arches cluster 

In the present section we confront the concept of F-dependent 
mass-loss rates with observations. To this end we use the large, 
and well studied sample of very massive stars in the Arches clus- 
ter, near the Galactic centre. The core of this young massive 
cluster contains 13 extremely luminous WNh stars, and a similar 
amount of bright early-type O, and Of stars. We adopt the stel- 
lar parameters L, X^, T*, and M, as obtained by iMartins et all 
J2008lfl in a comprehensive study of the Arches sample. 

In Sect. 15. II we start with a description of the properties of 
the Arches cluster stars. Using a fitting technique, we investi- 
gate in Sect. 15.21 whether the observed properties of the Arches 
stars are in line with a general F-dependent mass loss relation. In 
Se ct.l5.3lwe determine wind m asses from the mass-loss relations 
bv lGrafener & Hamannl(l2008l) . to test if the theoretical relations 
cover a realistic parameter range. 



5. 1 . Properties of the most massive stars in the Arches 
cluster 

The stellar parame ters of the Arches stars, as derived by 
IMartins et alj (I2008I) . are compiled in TablefTJ Spectral types 
reach from early O, and Of supergiants to late WNh subtypes. 
Stellar temperatures lie in the range = 30-40 kK, and lu- 
minosities in the range log(L*/L ) = 5.75-6.35, with the WNh 
stars showing systematically higher luminosities, and lower tem- 
peratures. Moreover, the large part of the WNh stars is H- 
deficient and N-enriched, with respect to solar values. Note, 
however, that some WNh stars display a solar hydrogen abun- 
dance at their surface. 

Based on the derived mass-loss rates, IMartins et al.l (I2008I) 
identify two distinct wind momentum - luminosity relations for 
the O/Of, and WNh stars. The WNh stars display systematically 
higher mass-loss rates than the O stars. This raises the question 
in which way the WNh stars differ from the O/Of stars. On the 
one hand, the different mass-loss properties could indicate a dif- 
ferent evolutionary stage, e.g., the WNh stars could be in the 
phase of core He-burning. On the other hand, parameters like 
surface abundances, or effective temperatures could be responsi- 
ble for the observed dichotomy. 



5.2. F-dependent mass-loss rates 

In this section we investigate whether the observed properties of 
the Arches cluster stars are in line with a F-dependent mass-loss 



4 N ote that the effective core temperatu re ?\, given by [Martins et aD 
d2008t) . is defined in the same way as bv lGrafener & Hamanrn2bo8l) . 
as the effective temperature related to the inner boundary radius R+ 
of the model atmosphere, i.e., by the relation L = 4nR 2 ^ crT*, where 
R+ is located at large optical depth. Due to the small density scale 
height in these layers, the precise value of the reference optical depth 
has alm ost no influenc e on the value of T+, so that the temperatures 
given bylMartins et all are ideally suited for the use with the relation by 
iGrafener & Hamannl . 
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Table 1. Stellar parameters of the most massive stars in the Arches Cluster. 



Star 


subtype 


T* 


L 


log(M) 




X H 


*N 




Mhmn 


M w2 


M w i 


pHeb 

e 






r w2 






[kK] 


[£o] 


r M0j 
yr J 






[%] 


[M a ] 


[Me] 


[M s ] 


[M e ] 










F8 


WN8-9 


33.7 


6.10 


-4.50 


1000 


0.18 


1.64 


37.1 


55.5 


51.0 


43.3 


0.63 


0.42 


0.46 


0.54 


F5 


WN8-9 


35.8 


5.95 


-4.64 


900 


0.22 


1.95 


29.5 


47.3 


36.5 


31.1 


0.58 


0.36 


0.46 


0.54 


F3 


WN8-9 


29.9 


6.10 


-4.60 


800 


0.27 


2.79 


37.1 


62.7 


63.2 


52.4 


0.68 


0.40 


0.40 


0.48 


F4 


WN7-8 


37.3 


6.30 


-4.35 


1400 


0.36 


2.10 


51.0 


93.6 


75.7 


66.4 


0.83 


0.45 


0.56 


0.64 


F2 


WN8-9 


34.5 


6.00 


-4.72 


1400 


0.40 


1.43 


31.8 


63.3 


49.2 


41.6 


0.68 


0.34 


0.44 


0.52 


F7 


WN8-9 


33.7 


6.30 


-4.60 


1300 


0.43 


1.86 


51.0 


101.2 


102.3 


86.3 


0.88 


0.44 


0.44 


0.52 


F6 


WN8-9 


34.7 


6.35 


-4.62 


1400 


0.54 


1.14 


55.4 


121.8 


119.2 


101.0 


0.97 


0.44 


0.45 


0.53 


F12 


WN7-8 


37.3 


6.20 


-4.75 


1500 


0.54 


2.26 


43.4 


97.5 


82.1 


70.0 


0.87 


0.39 


0.46 


0.54 


Bl 


WN8-9 


32.2 


5.95 


-5.00 


1600 


0.69 


2.41 


29.5 


82.0 


60.4 


49.8 


0.80 


0.29 


0.39 


0.47 


Fl 


WN8-9 


33.7 


6.30 


-4.70 


1400 


0.69 


1.45 


51.0 


133.0 


119.4 


100.9 


1.03 


0.40 


0.44 


0.52 


F9 


WN8-9 


36.8 


6.35 


-4.78 


1800 


0.69 


1.46 


55.4 


143.6 


131.3 


111.3 


1.07 


0.41 


0.45 


0.53 


F14 


WN8-9 


34.5 


6.00 


-5.00 


1400 


0.69 


0.49 


31.8 


87.4 


64.8 


54.0 


0.83 


0.30 


0.41 


0.49 


F16 


WN8-9 


32.4 


5.90 


-5.11 


1400 


0.69 


1.46 


27.4 


76.9 


56.0 


45.9 


0.76 


0.27 


0.37 


0.46 


F10 


04-6If 


32.4 


5.95 


-5.30 


1600 


0.69 


0.39 


29.5 


82.0 


69.1 


55.3 


0.80 


0.29 


0.34 


0.42 


F15 


04-6If 


35.8 


6.15 


-5.10 


2400 


0.69 


0.49 


40.1 


107.0 


96.9 


79.7 


0.93 


0.35 


0.38 


0.47 


F18 


04-61 


37.3 


6.05 


-5.35 


2150 


0.69 


0.39 


34.3 


93.3 


82.5 


66.9 


0.86 


0.32 


0.36 


0.44 


F20 


04-61 


38.4 


5.90 


-5.42 


2850 


0.69 


0.30 


27.4 


76.9 


57.4 


46.7 


0.76 


0.27 


0.36 


0.45 


F21 


04-61 


35.8 


5.95 


-5.49 


2200 


0.69 


0.39 


29.5 


82.0 


70.2 


56.0 


0.80 


0.29 


0.33 


0.42 


F22 


04-61 


35.8 


5.80 


-5.70 


1900 


0.69 


0.39 


23.7 


68.0 


52.5 


41.3 


0.70 


0.24 


0.32 


0.40 


F23 


04-61 


35.8 


5.80 


-5.65 


1900 


0.69 


0.69 


23.7 


68.0 


51.5 


40.7 


0.70 


0.24 


0.32 


0.41 


F26 


04-61 


39.8 


5.85 


-5.73 


2600 


0.69 


0.40 


25.5 


72.3 


56.6 


45.0 


0.73 


0.26 


0.33 


0.42 


F28 


04-61 


39.8 


5.95 


-5.70 


2750 


0.69 


0.40 


29.5 


82.0 


71.5 


56.8 


0.80 


0.29 


0.33 


0.41 


F29 


04-61 


35.7 


5.75 


-5.60 


2900 


0.69 


0.30 


22.1 


64.1 


44.8 


35.6 


0.67 


0.23 


0.33 


0.42 


F32 


04-61 


40.8 


5.85 


-5.90 


2400 


0.69 


0.29 


25.5 


72.3 


59.3 


46.6 


0.73 


0.26 


0.31 


0.40 


F33 


04-61 


39.8 


5.85 


-5.73 


2600 


0.69 


0.39 


25.5 


72.3 


56.6 


45.0 


0.73 


0.26 


0.33 


0.42 


F34 


04-61 


38.1 


5.75 


-5.77 


1750 


0.69 


0.40 


22.1 


64.1 


46.0 


36.3 


0.67 


0.23 


0.32 


0.40 


F35 


04-61 


33.8 


5.70 


-5.76 


2150 


0.69 


0.20 


20.6 


60.5 


43.1 


33.6 


0.64 


0.22 


0.31 


0.39 


F40 


Q4-6I 


39.8 


5.75 


-5.75 


2450 


0.69 


0.40 


22.1 


64.1 


44.5 


35.5 


0.67 


0.23 


0.33 


0.42 



Notes. Designations, subtypes, stellar te mperatures (T+), lumi nosities (L), mass-loss rates (M), terminal wind velocities ([),*,), hydrogen, and nitro- 
gen surface mass fractions (X S H , X* ), from lMartins et al.ld2008h . and mass estimates from the present work: He-burning mass es (MH e b), and homo- 
geneo us masses (M hom ) according to Eqs. ll llfTJl with coefficients from rows 11, & 16 in Table lA. fl wind masses according to l Grafener & Hamannl 
(2008) for Z = 2Z Q (M W 2), and Z = 1 Z Q (M w i). The corresponding Eddington factors T c are given in the last four columns. 

sible to constrain the free parameters Mo, fr, fr, and fj by a 
A^-fit, based on the observed values of M, L, T*, and X^. 

The hydrogen surface abundance plays a crucial role in 
this fitting process, as it enters the estimate of F e via Xe (Eq.[8]l, 
and via Mh om (L, Xh). The wide spread in surface abundances for 
the WN stars in our sample is important for the verification of an 
actual T-dependence. Without variations in X^, a F-dependence 
could not be easily distinguished from an L-dependence, if there 
exists an underlying M(L) relation that does not depend on other 
parameters. In the following we will see that this may indeed 
happen for the O stars within our sample, which all have the 
same X^. 

Using Eq.[l4] we perform ^ 2 -fits to the observed mass-loss 
rates of the sample stars, based on the observed values of L, T*, 
and X^. Dependent on whether T e is estimated on the basis of 
Mueb(L), or Mhom(L, Xh), we designate the obtained values of fr 
in Eq.[14]as f^ om , or /j? eb . To extract the physically relevant pa- 
rameters we perform a series of tests where the fitting parameters 
fr, fh, and fj are partly set to zero. 

We start our investigation using only the 13 WNh stars in 
the Arches sample. The results of our ;f 2 -fits are compiled in 
Table|2] We start with the most simple case of a constant mass- 
loss rate, i.e., with Mo as the only free parameter. The obtained 
X 1 = 0.544 obviously marks the lower end of the achievable 
fit quality. In the next four rows of Table|2] we allow for one 
more free parameter apart from Mo- Among those four cases the 
best result is obtained for a T-dependent mass-loss relation under 



Table 2. Results of various^ 2 fits of Eq.[14]to the observed mass- 
loss rates of the 13 WNh stars in the Arches sample. 



X 1 


M 


f horn 

Jr 


fHeb 

Jr 


ft 


fr 


0.544 


-4.721 - 


0.108 


-3.617 


2.578 








0.536 


-4.752 




-0.330 






0.368 


-4.820 






0.720 




0.531 


-4.751 








1.204 


0.354 


-4.800 






0.843 


-1.438 


0.039 


-2.802 


4.208 




-0.845 




0.098 


-3.537 


2.703 






-1.056 


0.037 


-5.244 




-2.937 


1.749 




0.496 


-4.859 




-0.787 




2.467 


0.037 


-2.768 


4.293 




-0.920 


0.495 


0.035 


-5.262 




-3.008 


1.724 


0.577 



Notes. Parameters indicated by (-) have been set to zero. 

relation. For this purpose we adopt a general mass loss relation 
in the form of a power law with four free parameters 

log(M/ M Q yr 1 ) = M + fr ■ log(F e ) 

+ Sl ■ (logCL/L ) - 6.0) + f T ■ (log(7VK) - 4.5) . 

To estimate T e , we use Eq.[8] with the relations for 
Mhom(L, Xu), and Mn e \,{L) from Sect. 14.11 This way, it is pos- 
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Fig. 3. Under the assumption of chemical homogeneity, the 
Arches cluster stars display a pronounced F-dependence. 
Squares indicate stars with spectral subtypes WNh (filled blue: 
H-deficient; unfilled: normal H abundance), and circles O/Of. 
The black dashed line indicates the fitted mass-loss relation 
CEq.n~4b for WNh stars (parameters from Table[2]row 2). The 
red dotted line indicates the corresponding fit for the complete 
sample (parameters from Table[3]row 2). 
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Fig. 4. Under the assumption of central He-burning, the Arches 
cluster stars display no general T-dependence. Symbols are the 
same as in Fig. [3] the black dashed line indicates a fitted relation 
(Eq.fPfl) for WNh stars (parameters from Table|2]row 3). 



the assumption of homogeneity. Note the big difference in x 1 
between the homogeneous assumption, and the assumption of 
core He-burning. This is illustrated in Figs. [3] & |4] where the 
sample stars indeed form a well pronounced mass-loss relation 
under the assumption of homogeneity (Fig.|3]l, but not under the 
assumption of He-burning (Fig.|4j». At this point we note that the 
two plots represent two extremes, and that the real M-T relation 
may lie in between the two cases. 

The qualitative difference between the two plots originates 
from the fact that the adopted M-L relation in Fig. [3] depends 
on X^, while the relation in Fig. [4] does not. In both plots the 
O stars form a separate sequence that merges into the WR se- 
quence. As explained above, this does not necessarily indicate 
a true F-dependence for the O star sample, as these stars show 
no variations in X^. The WNh stars, on the other hand, display a 
substantial spread in X^. Notably, in Fig. [4] the H-deficient stars 
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Fig. 5. Fitted mass-loss relation for the complete Arches sample, 
under the assumption of chemical homogeneity (Eq.[l5j red dot- 
ted line). Symbols are the same as in Fig. [3] The observed mass- 
loss rates are scaled to a temperature of log(r*/K) = 4.5, based 
on the temperature dependence in Eq.[l5] w ith M oc t~ 4116 . For 
comparison, we plot the relation bv lGrafener & Hamannl (120081) 
for the same temperature, log(L) = 6.3, Xu = 0.7, and Z = Z 
(black solid), as well as Z = 2Z (grey solid). 

Table 3. Results of various^- 2 fits of Eq.fT4lto the observed mass- 
loss rates of the complete Arches sample. 





x 2 


Mo 


j~hom 

Jr 


ft 


h 




6.557 


-5.190 










0.942 


-2.871 


4.550 








2.037 


-5.179 




2.068 






4.160 


-4.880 






-8.458 




0.783 


-1.694 


6.872 


-1.239 






0.457 


-2.957 


3.997 




-4.116 




0.429 


-2.420 


5.091 


-0.557 


-3.751 



Notes. Parameters indicated by (-) have been set to zero. 



show a significant shift towards lower F e , i.e., a 'hook' in the 
M-T relation. This can be explained by the dependence of T e on 
X^ in Eq.[8] with T e oc 1 + X^. The fact that this effect is pre- 
cisely compensated in Fig. [3] suggests that the internal structure 
of the WNh stars in our sample is close to homogeneity. More 
precisely, it shows that there exists a relation between core, and 
surface hydrogen abundance. A similar effect may thus occur for 
inhomogeneous stars that have a similar degree of inhomogene- 
ity, so that the differential behaviour of the homogeneous case 
is preserved. However, also in this case, the largest part of the 
sample stars need to be in phase of core H-burning to display 
the observed behaviour. Only few individual objects, with the 
highest He surface enrichment, may already have reached the 
He-burning phase. 

In the remainder of Table|2] we allow for three or four free 
parameters. For these cases similar fitting results are obtained 
for the He-burning, and the homogeneous case. However, the 
obtained values for /j? eb are all negative. Such a relation would 
imply that mass loss ceases close to the Eddington limit, which 
seems to be rather unphysical. An inspection of Fig. [4] shows 
that the assumption of He-burning would indeed imply a slightly 
negative slope of the mass-loss relation for WNh stars. The fact 
that for this case a similar fit quality can be achieved as for the 
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homogeneous case is presumably a sign that the number of fit- 
ting parameters becomes too large for the small sample. 

In Fig.[3]it is notable that our fitted M(T e ) relation for WNh 
stars (black dashed line in Fig. [3} seems to be slightly detached 
from the rest of the sample, similar to the two d ifferent mass- 
loss re lations for WN and O/Of stars identified bv lMartins et alJ 
d2008l) . Fig.El however, suggests that there might exi st a smooth 
transi tion between the WNh, and O/Of stars (cf. IVink et all 
1201 ll) . To investigate this possibility we have performed the 
same fitting procedure as previously, for the complete Arches 
sample. Because the O stars are unlikely to be in the phase of 
core He-burning, we have however only focused on the homo- 
geneous case. 

The results are combined in Table[3] Again, we start with the 
simplest case of a constant mass-loss rate Mo with^- 2 = 6.557. 
The next three rows allow for one more free parameter. Again, 
a F-dependent mass-loss relation is strongly favored with x 2 = 
0.942. Interestingly, the following fits with more free parameters 
strongly favor a temperature dependence wi th M oc T~ 4 , which is 
very si milar to the dependence predicted by Gra fener & Hamannl 
(120081) . The two fits with the lowest^- 2 in Table[3]do generally 
support the key features of their mass-loss predictions, which are 
1) a strong F-dependence, 2) a weak luminosity dependence, and 
3) a strong increase in mass loss for decreasing T*. 

From row 6 in Table[3]we obtain the mass-loss relation 



logtM/Moyr- 1 ) = 

- 2.957 + 3.997 ■ log(r e ) - 4.116 ■ (log(7VK) - 4.5) . 



(15) 



In Fig.|5]we compare this relation with observed mass-loss rates, 
compensating for the temperature dependence. Notably, the pre- 
vious distinction between the O/Of, and the WNh sample from 
Fig.|3]does not exist anymore in this plot. The temperature differ- 
ence between these two samples can thus account for the differ- 
ences in their mass-loss rates. We note however that the depen- 
dence on T+ is uncertain, as it only follows from the combined O 
+ WR sample. From the WR sample alone there is no evidence 
for such a dependence, presumably due to the small spread in 
T+ . At this time it is not clear at which point the transition to a 
T-dependent mass loss relation occurs, and to which extent such 
a relation is applicable to (parts of) the O star sample. 

5.3. Wind masses for the Arches cluster stars 

In the previous section we have shown that the observed prop- 
erties of the Arches cluster stars are in line with a T-dependent 
mass-loss relation. The comparison in Fig. [5] shows that, under 
the assumption of chemical homogeneity, the resulting mass- 
loss relation is in good ov erall agreement with the relation by 
iGriifener & Hamannl (120081) for a high metallicity of Z ~ 2 Z . In 
the prese nt section we use win d masses that follow from the re- 
lation by Gra fener & Hamannl to investigate if the Arches stars 
conform with this mass-loss relation, and if the obtained wind 
masses cover a plausible range. 

We derive wind masses by computing the Eddington factor 
r e that is needed to explain the observed mass-loss rates M, for 
given stellar parameters L, T*, Xh, and Z, from Eqs.[3}{5] The 
wind masses follow from T e , according to Eq.[8] In TableQ]we 
list the resulting masses, and Eddington factors. The different 
columns in TableQ]denote values obtained for an adopted metal- 
licity of 2Z (M W 2), and 1Z (M w i). The results are compared to 
homogeneous masses, and He-burning ma sses from S ect. |4.1l 

We note that the mass-loss relation by iGrafener & Hamannl 
(120081) is obtained for stellar winds with an optical depth of the 
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Fig. 6. Mass-luminosity relations L(M) for the most massive 
stars in the Arches cluster, for an adopted metallicity of Z = 
2Z . Plotted are observed luminosities, vs. different mass es- 
timates. Red diamonds indicate L(Mh om ), with homogeneous 
masses Mh om (L,Xn) computed from Eqs.Qj] &[T2] with coef- 
ficients from row 1 1 in Table lA. ll Black/blue symbols indicate 
L(M W i n d), with wind masses obtained from the mass-loss pre- 
scription bylGrifener&HamanB (Hool) (Eqs.|3}S forZ = 2z o- 
For the latter, we distinguish between WNh subtypes (squares), 
and O/Of subtypes (circles). Among these, blue filled symbols 
indicate H-deficient (WNh) stars. The dashed lines indicate the 
He-MS (left), and the ZAMS (right) according to our relations 
in Sectltll 
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Fig. 7. Mass-luminosity relations L(M) for the most massive 
stars in the Arches cluster, for an adopted metallicity of Z = 1Z Q . 
Symbols are the same as in Fig. [6] 



order of one, or higher, and that the O stars with lower luminosi- 
ties in our sample may just lie outside the applicable range of 
this relation. Nevertheless we include them here, also to check 
the validity of the mass-loss relation. 

In Fig. [6] we compare the resulting M-L relation for Z = 
2Z Q , with the expected relation for homogeneous stars. We 
find a notable overall agreement between both M-L relations. 
Morphologically, both relations show two branches in the M— 
L diagram. The first branch, with higher luminosities for given 
wind mass, consists of H-deficient stars with subtypes WN7-9h 
(blue filled symbols). The second branch, with lower luminosi- 
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ties, consists of stars with normal hydrogen surface composition 
(empty black symbols). This group covers spectral subtypes 04- 
61, 04-6If+, and WN8-9h. Our theoretically predicted masses 
are indicated by red diamonds. Also they show two branches, 
one on the ZAMS (at lower luminosities), and one in the region 
of our previous high-luminosity branch. 

For the theoretically predicted masses Mh om (L, X^) (red dia- 
monds) the interpretation is very simple. By definition, homoge- 
neous stars with solar H-abundance populate the ZAMS, while 
H-deficient stars have a higher mean molecular weight, and thus 
display higher luminosities. 

For the wind masses M W 2 (black/blue symbols), the differ- 
ence between both branches mainly originates from the higher 
mass-loss rates of the WNh stars, and the temperature differ- 
ence between the WNh, and the O stars. For the H-deficient 
branch, which consists only of WNh subtypes, the derived wind 
masses are in very good overall agreement with the homoge- 
neous masses. For the H-rich branch, the derived wind masses do 
not fit the homogeneous masses, but are systematically lower, as 
expected for stars that have just evolved away from the ZAMS. 
Again, we have to keep in mind that the stars with the lowest 
luminosities in this group, may lie beyond the applicable regime 
of the theoretical mass-loss relation. 

For the H-deficient WNh stars, we only find a good quanti- 
tative agreement between empirically determined wind masses, 
and theoretically predicted homogeneous masses, if we assume 
a high metallicity of 2Z . A similar metallicity is indeed mea- 
sured indirectly for the WNh stars in the Arches sample, based 
on their nitrogen abundances. However, more recent investiga- 
tions of stars in the GC, favor a solar value (cf . Sect. 1 6. U . Using 
the mass-loss relation bv lGrafener & Hamannl (120 08) with solar 
Z, the derived masses are lower, but the morphological similar- 
ities remain (cf. Pig-EJ. In this case the wind masses of both 
branches in the diagram imply an inhomogeneous stellar struc- 
ture. 

We conclude that, usi ng the mass-loss prescription by 
Gra fener & Haman nl (120081) . we find two branches in the M-L 
diagram. The more luminous branch is populated by H-deficient 
WNh stars. Based on an adopted metallicity of 2Z Q , the derived 
wind masses for these stars would be in agreement with chemi- 
cally homogeneous stars. Interestingly, these stars also show the 
strongest nitrogen enrichment (cf. Table[TJ, which points to a 
previous evolution with episodes of strong mixing, and/or strong 
mass loss, that could indeed lead to a quasi-homogeneous struc- 
ture of these objects (cf. Sect. l6.4l i. The H-rich stars, on the other 
hand, populate a branch with lower luminosities, in agreement 
with stars that have undergone a normal evolution off the ZAMS. 

At this point we want to note that our p reference for a high Z 
could mean that the mass-loss relation by iGrafener & Hamannl 
J2008h still underestimates the true mass loss rates. In this case 
the high Z would compensate for lacking elements in the mod- 
els. Another possibility would be, that the sample stars are not 
homogeneous, and thus F e is underestimated by our approach. 
However, in any case the wind masses cover a plausible parame- 
ter range. Notably, almost all of the sample stars lie below the He 
main-sequence in the M-L diagram. In agreement with our find- 
ings from Sect. 15. 21 the largest part of the Arches cluster stars is 
thus most likely still in the phase of core hydrogen burning. 

6. Discussion 

In the present work we have investigated the mass-loss proper- 
ties of the most massive stars in the Arches cluster. We have used 
a semi-empirical method that is based on observations, or more 
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Fig. 8. Mass-luminosity relations L(M) for the most massive 
stars in the Arches cluster, for an adopted metallicity of Z = 2Z Q . 
To investigate the influence of systematic uncertainties on our re- 
sults, the observed luminosities have been artificially increased 
by 0.15 dex. To derive wind masses, the mass-loss rates have 
been increased according to M oc L 3 / 4 , as expected for the results 
from recombination line analyses (cf., lHamann & Koesterkd 
119981) . The symbols are the same as in Fig.l6l 



precisely, on the results of spectral analyses, and on theoretically 
predicted M-L relations. As a result we found evidence for the 
existence of a largely F-dependent mass-loss relation, and ob- 
tained information about the evolutionary status of the Arches 
cluster stars. 

In the following we discuss our results. We start with the 
question how the results are affected by potential uncertainties 
in the observed stellar parameters (Sect. 16. it . In Sect. 16. 21 we 
discuss the implications for the mass loss properties of massive 
stars, and in Sect. [63] we reconsider the potential to obtain pre- 
cise wind masses. Finally, in Sect. 16.41 we discuss the implica- 
tions on the evolutionary status of the Arches cluster stars. 

6.1. Observational uncertainties 

Here we discuss to which extent our results are affected by un- 
certainties in the stellar parameters, that we use as input for our 
method. We give an overview of potential error sources, discuss 
their quantitative importance, and their influence on the qualita- 
tive outcome of our work. Generally it can be said that, due to 
the steepness of the proposed M-Y relation, our results are rel- 
atively insensitive to such uncertainties. In addition, most error 
sources are of systematic nature, and hardly affect the qualitative 
out come of our work. 

iMartins et al] d2008l) quote uncertainties of +3 kK (2 kK) on 
T eS for O stars (WNLh), +0.2 on log(L), 0.2 (0.1) on log(M), 
and ±50% (+30%) on abundances. If these (rather conservative) 
error estimates would be of purely statistical nature, they would 
be devastating for our results. The fact that we actually find well- 
defined relations for the sample stars, indicates that potential er- 
rors are chiefly of systematic nature. The possibility that there 
are differences in the systematics between O stars, and WNh 
stars, could however influence our results, particularly regarding 
the tempe rature dependen c e in Se ct. 15.21 

Also | Crowfher et al.l (1201 Oh recently pointed out that 
IMartins et al.l (12008 ) may have systematically under-estimated 
the stellar luminosities of the Arches cluster stars, by using 
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a relatively small dis tance estimate (m - M = 14.4 + 0.1, 
lEisenhauer et al.ll2.005h. an d a low foreground extinction (Ak = 
2.8 + Q.llStolteetalj2002l) for t he GC. Using m-M = 14.5 ±0.1 
dReidll993h .andAi y = 3.1+0.1 dKim et alJ2006hlCrowther eTaT] 
find slightly higher luminosities for the two brightest stars in the 
sample (log(L/L Q ) = 6.5 instead of 6.35). 

Although a systematic increase of log(L) by 0.15 dex would 
be significant, the T e , as derived from the mass-loss relations 
in Sect. [3] would hardly be affected. The reason is that spectro- 
scop ic mass-loss rates der ived from recombination line analyses, 
as in lMartins et al.l(l2008l) . scale with M oc L° 75 . As the luminos- 
ity scaling in our mass-loss relations is very similar (Eqs.|6]|7] & 
Eq.0, the derived T e would nearly stay the same. This would 
however still result in uncertainties in the derived wind masses, 
basically with M w ; n[ j oc L(l + X^) (cf. Eq.0. The effect of a sys- 
tematic luminosity increase by 0.15 dex is illustrated in Fig. [8] 

For the mass-loss relations in Sect. 15. 21 a luminosity increase 
mainly results in a shift of the inferred mass-loss rates with M oc 
L 75 . In addition, the change of the slope of the M-L relation 
changes the slope of the mass-loss relations, with respect to T e . 
E.g., based on relation No. 1 in Table lA. fl an increases of 0.15 
in log(L) for a star with log(L/L ) = 6.3, and X = 0.7 (i.e., 
log(r e ) = -0.386) causes an increase of T e by 0.048 dex. For a 
star with log(L/L ) = 5.7 (log(r e ) = -0.647), the same increase 
causes F e to increase by 0.073 dex. Consequently, the exponents 
in our empirical mass-loss relations increase by a factor of 1 + 
(0.073 -0.048)/(0.647- 0.386) = 1.1. 

Additional uncertainties in the mass-loss rates are due to 
wind clumping (with M oc JJa), and the unknown wind ve- 
locity field. Again, these would introduce a shift in M, but only 
moderately affect the derived T e , because of the steepness of 
the proposed T-dependence. Based on the aforementioned un- 
certainties, we thus expect systematic effects on our quantita- 
tive results, particularly concerning mass-loss rates, and wind 
masses, but only moderate relative changes, i.e., our results stay 
qualitatively the same. 

Apart from that, the metallicity Z is of major importance. In 
contrast to the above, Z does affect the F-dependence in our the- 
oretical mass-loss relations, but not the observed properties of 
the sample stars. As discussed in Sect. 15.31 it thus affects the de- 
rived wind masses. In the relation bv lGrafener & Hamannl d2008l 
Eq.FJ), the metallicity determines the value Fo(Z) (Eq.|4j», which 
roughly marks the value of T e for which the mass loss starts to 
increase due to the proximity to the Eddington limit. Moreover, 
the e xponent fi(Z) defines the steepness of the mass-loss relation. 
Also lVink & de Kotej d2005l) find a mass-loss dependence with 
M oc Z 0M for a typical Wolf-Rayet model. The reason for the 
strong Z-dependence lies in the dominant role of Fe line opaci- 
ties in the radiative acceleration of hot star winds. The Z depen- 
dence in the wind models thus chiefly reflects a dependence on 
the Fe abundance. 

As mentioned before, the metallicity of the Arches clus- 
ter is still uncertain. It has been determined indirectly, on 
the b asis of observed nitrogen surface abun dances of WNh 
stars dNajarroet al.ll2004t IMartins et aT]l2008l) . These stars are 
strongly nitrogen enriched, obviously by material which has 
been processed in the CNO cycle in the stellar core. This ma- 
terial is exposed at the stellar surface, presumably by mixing 
and/or mass loss. The observed N abundance thus provides a 
lower limit for the initial abundance of C+N+ O. For the WNh 
star s in the Arches cluste r iNaiarro et al.l (120041) find X N « 1 .6%, 
and IMartins et alJ d2008l) find an average of X^ as 1.7%, with 
values reaching up to X^ = 2.79%. Both authors conclude that 
the initial metallicity of the Arches cluster is only slightly super- 



solar. However, taking into account the recent d ownward revi- 
sion of the s o lar oxygen abundan ce (according to lAsplund et al.l 
I2005L 120091: iPereira et alJ l2009l it follows that X C+N+0 , Q * 
0.8%), a metallicity of 2Z seems more likely. As the solar Fe 
abundance has hardly been affected by this revision, a simple 
scaling would imply that Xp e /Xp e<Q = 2. 

We note howev er that recent spect ral analyses of LB Vs 
dNajarro et al ,l2009h . and red supergiants dDavies et alj200Sj). as 
well as earlier studies of cool stars in the GC region dCarr et all 
l2000HRamrrez et al.l2000l:ICunha et alj2007l) indicate a solar Fe 
abundance, and an increased abundance of alpha elements. The 
resulting uncertainty in the derived wind masses is of the order of 
20% (see TableQ]). Generally, radiatively driven winds become 
stronger for higher Z, i.e., the derived wind masses become sys- 
tematically larger if a higher Z is adopted. As the derived mass 
ratios are hardly affected, we do not expect a qualitative impact 
on our results from Sect. 15.31 The preference for a high Fe abun- 
dance in this work, particularly in the comparison in Fig. [5] may 
indicate that either the T e in Sects. [5721 and 15.31 are underesti- 
mated, due to the assumption of homogen eity, or the mass-loss 
predictions by IGrafener & Hamannl d2008l) are still too low. In 
any case, we do not expect a qualitative impact on our results, 
as the differential stellar properties within the sample are hardly 
affected. 



6.2. The mass-loss properties of very massive stars 

The main result of the present study is the empirical con- 
firmation of a T-dependent mass-loss relation for stars ap- 
proaching the Eddingt on limit. This result is consistent with 
theoretical prediction s ( Grafener & Hamannl 120061: IVinkl 120061 
IGrafener & Ham ann 2008b IVink et al.ll201 lh ! as summarized in 
Sect.|3] A similar F-depe ndent relation has also been proposed 
for the winds of LBVs, bv lVink & de Ko ter (2002). When we ex- 
press our empirical results in the form M(F e , T*, L, X^), we find 
that the dependencies on L, and X^ are very weak, indicating 
that r e , and potentially T*, are the physically most relevant pa- 
rameters. The te mperature dependence in Eq . fTSl is in remarkable 
agreement with IGrafener & Hamannl (120081) . however, based on 
the present data, this result seems not fully established. 

The empirical confirmation of a T-dependence marks a 
paradigm shift in the way we think about the development of 
stellar mass loss through subsequent evolutionary phases. The 
physical effect that causes (LBV and) WR-type mass loss is the 
proximity to the Eddington limit. During their evolution, the 
cores of massive stars become chemically enriched, and their 
L/M ratio increases. When Eddington factors of order unity are 
reached, a strong WR-type stellar wind will develop naturally. 
With the removal of the outer envelope, in combination with 
mixing processes, a hydrogen deficiency will develop at the stel- 
lar surface. It is therefore not the hydrogen deficiency that de- 
fines the onset of the WR phase, as is usually assumed in stellar 
evolution calculations. Rather, the WR-phase starts as soon as 
Eddington factors of order unity are reached. 

This effect will have strong impact on the evolution of mas- 
sive stars, in terms of conventional mass loss, but also in terms 
of loss of angular momentum (e.g.. lLangerlll998l) . For fast ro- 
tating stars, the Eddington parameter effectively increases near 
the equator, due to centrifugal forces. The eff ect of rotation on 



mass l o ss has been studied for OB stars, e.g. by Friend & Abbott 
(|1986|) : [Pelupessv et al I d2000l) : ICure & Riall d2004h : ICure et al. 
d2005h : Madura et all d2007h . Close to rotational breakup (i.e., 
close to the £2-limit), the latter two of these works find shallow 
wind solutions, for which the mass loss does not significantly ex- 
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ceed the mass loss w ithout rotation (only by a factor of two, see 
IMadura et al.ll2007h . For WR-type winds, the effect of rotation 
has not yet been investigated. The strong F-dependence however 
implies that the combined effect of rotation and high F may still 
be efficient at high rotational speed (i.e., close to the FQ-limit). 
This effect may be of paramount importa nce for fast rotating WR 
stars, like the progenitors of long GRBs (lYoon & Langer 20051) . 
At the present stage this is however highly speculative. 

In addition to that, the continuous transition between the 
O/Of, and WNh stages (cf. Fig.[5]l has consequences for our 
theoretical understanding of the winds of very massive stars. 
Notably, this suggests that they have the same driving mecha- 
nism, namely radiation pressure. Nevertheless, there are impor- 
tant differences between O-star winds and WR-winds. While the 
former are expected to be optically thin at the sonic point, the 
latter are optically thick, i.e., the wind is initiated and already 
accelerated in part below the surface. The transition between 
these tw o regimes has be en resolved by the model computa- 
tions of IVink et al.l (1201 lh with their different r-dependen cies 
for the OB, and the WR regime (Eqs.|6] & |7). IVink et al.l find 
that this transition occurs for models with a wind efficiency 
77 = MvcoliL/c) « 1, which roughly corresponds to the point 
where the winds become optically thick. 

The properties of optically thick, WR-type winds have been 
suggested to be determined by the so-called Fe-opacity peaks 
in deep atmospheric layers, close to the sonic point (e.g. , 



Pistinner & Eichlerll995l:lHeger & Langerfl996USchaeredl996t 



Nugis & Lamersll2002l). This has been conf irmed by WR wind 
models of lGrafener & Hamannl (120051 12008). who proposed that 
the formation of optically thick winds is supported close to the 
Eddington limit, because the density scale height in the deep at- 
mospheric layers increases. In contrast to this, OB star winds are 
thought to be dominated by the outer wind physics. Within th e 
classical theory of radiatively driven winds (ICastor et al.ll 1975k 
their wind properties are determined at a critical point with a 
much higher wind speed, which is typically of the order of the 
escape speed. The wind energy at this point is thus of the same 
order of magnitude as the gravitational energy. Small changes 
in the (effective) gravitational potential thus only have small in- 
fluence on the mass loss of OB stars, i.e., the F-dependence is 
expected to be weak. The fact that we find observational evi- 
dence for a strong T-dependence of WR mass loss, in agreement 
with the model predictions, shows that there is indeed a substan- 
tial difference in the way the mass-loss rates of OB and WR stars 
are determined. 

A strong temperature dependence of WR-type mass loss 
would support the i mportance of the deep layer s for the wind 
physics of WR stars. iGrafener & Hamarmn feOOS) predict a sim- 
ilar dependence, with M oc T 1 ^ 3 . They explain this depen- 
dence by the necessity that the critical point of the equation of 
motion needs to be located in a specific temperature regime, 
and thus at higher optical depths for cooler stars. They also 
show that this temperature dependence is in agreement with 
the observed line strengths of galactic WNL stars (cf. Fig. 2 in 
IGrafener & Hamanr]l2008l) . Quantitatively, the temperature de- 
pendence might be of major importance because it leads into 
phases of extremely high mass loss, potentially even to a smooth 
transition into the LBV phase. We note, however, that our em- 
pirical result arises from a comparison of the hotter O/Of stars 
vs. the cooler WNh stars. Particularly for the low-luminosity end 
of the O star sample we concluded in Sect. l5.3l that they might 
have a systematically different internal structure than the WNh 
stars. Moreover, for the relatively weak winds of these objects, 
the concept of T-dependent mass loss might be questionable. The 



fact that we actually find a relation that includes these objects is 
thus rather surprising. Because of these concerns it would be de- 
sirable to back up this result by complementary studies of WNh 
stars in different temperature regimes. 

In Fig. [5] we compare our emp irical relation (Eq.[T5jl with 
theoretical mass-loss predictions bv lGrafener & Hamannl d2008h 
for different metallicities Z. For a high Z of 2Z , we find a 
good quantitative agreement between theory and observation. 
Qualitatively, the slope of the empirical F-dependence is very 
we ll reproduc e d by t he models. The same holds for the relation 
bv lVink et al.1 (1201 ll) . Their exponent of 3.99 in Eq. [7] matches 
the exponent in Eq. [15] very precisely. The major difference be- 
tween the model predictions in Sect.[3]thus lies in the strong tem- 
peratu re dependence that is predicted by IGrafener & Hamannl 
d2008l 

Summing up our results, we expect that massive stars close 
to the Eddington limit tend to form WR-type winds with distinct 
properties from classical radiatively driven winds. In the WR 
regime we find mass-loss relations of the form M(F e , T+, L, X^) 
which fundamentally depend on the Eddington parameter T e . 
Our theoretical models predict 1) a strong F-dependence, 2) 
a weak dependence on L, and X ^, and 3) a strong t emper - 
ature dependence, according to Gr afener & Hamannl d2008l) . 
Empirically, point 1), and 2) are convincingly confirmed in our 
present study. The important temperature dependence 3), is in 
agreement with our present results, however, demands for fur- 
ther empirical confirmation. 

6.3. Wind masses for very massive stars 
In Sect.15.31 we used the mass-loss relations by 



IGrafener & Hamannl (|2()08) to derive the present masses 
of the Arches cluster stars, based on their observed mass-loss 
rates. The potential to estimate stellar masses in this way is 
a very interesting aspect of this work, in particular because 
our mass estimates are fundamentally different from wind 
masses that have pre viously been obtained for OB stars (e.g., 
iKudritzki et aDfT992h . 

In Sect. 16. 21 we already discussed that the dynamics of clas- 
sical radiation-driven winds is dominated by the outer part 
of the wind, while WR-type winds are strongly influenced 
by deep atmospheric layers. Classical OB stars thus show a 
well-defined relation between their escape velocity at the stel- 
lar surface, and their t erminal wind velocity, t hat helps to de- 
termine their masse s ( Kudritzki & Pulsl l2000t iPrinia & Massal 
119981: lLamers et al.lll995l) . Only IVink & de Koted d2002h previ- 
ously used a F-dependent mass-loss relation to constrain the stel- 
lar mass of the LBV AG Car, however, their models still depend 
on the outer wind physics, and thus on the ratio i>oo/u e sc- 

WR-type winds depend mostly on the inner wind physics. 
We have seen that this leads to a steep T-dependence for stars 
that are close to the Eddington limit, i.e., we are capable to de- 
termine precise F e , without consideration of the terminal wind 
velocity. This method is very promising, because the masses of 
stars close to the Eddington limit are already well known, sim- 
ply because F e is of the order of one. The wind models thus 
only need to provide an 'improved guess' of r e . Moreover, ob- 
servational uncertainties are of minor importance because of 
the steepness of the F-dependence. The main uncertainties arise 
from the models, chiefly from the contribution of metal lines 
to the 'physical' Eddington factor T(r). However, even if this 
contribution would be systematically wrong, the models still 
provide reliable mass ratios, i.e., the mass-loss relation can be 
calibrated with observations. In the present work we have per- 
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formed a first step towards such a calibration. It seems, however, 
that the remaining uncertainties in the metallicity, and the dis- 
tance/extinction towards the Arches cluster demand for comple- 
mentary studies in different environments. 

6.4. Evolutionary status of the Arches cluster stars 

As a by-product of our present study we obtain information 
about the internal structure, and thus about the evolutionary sta- 
tus of the WNh stars in the Arches cluster. Most importantly, 
we find direct evidence that the largest part of these objects is 
still in the phase of core hydrogen burning. In Fig. [4] we have 
shown that, if the stars were in the phase of core He-burning, 9 
of the 13 WNh stars in the sample would have Eddington factors 
T e > 0.75 (cf. also TableQ}. Taking into account that F(r) > T e 
(Eq.|2]l, this would imply that these stars are extremely close, 
or even above the physical Eddington limit. The reason for this 
are the high hydrogen surface abundances of these stars (up to 
solar). 

The M-L relations used in this work are based on the as- 
sumption of a simple, chemically homogeneous internal struc- 
ture of the sample stars. In Sects. lBT2l & l5.3l we find indications 
that this assumption may indeed be fulfilled by most of the WNh 
stars in our sample. However, due to uncertainties in the obser- 
vational data, and the models, we cannot draw firm conclusions 
at this stage. 

E.g., we found in Sect. 15. 31 that the wind masses of the H- 
deficient WNh stars are in good overall agreement with the 
expected masses of chemically-homogeneous stars (cf. Fig . [6j> . 
This result however depends on the adopted metallicity in our 
mass-loss relation. In Sect. 15.21 we found evidence for chemical 
homogeneity, based on mass ratios, which are much more certain 
(cf. Sects. |6*TI & 16.3b . Figs.[3]&|4]show that, under the assump- 
tion of chemical homogeneity, the relative masses of the sample 
stars adjust in a way that we observe a smooth M(F e ) relation. 
Under the assumption of core He-burning, i.e., the completely 
in-homogeneous case, we obtain differences in M by a factor of 
~10 for similar T e , particularly for H-deficient WNh stars. This 
suggests that most of the H-deficient WNh stars in the Arches 
cluster are close to homogeneity, although a similar behaviour 
may also originate from a stellar sample where all stars have a 
similar degree of inhomogeneity. 

A chemically homogeneous structure is expected for single, 
fast rot ating ma s sive st a rs, e.g., according t o evolutiona r y mod - 
els by, iMaederi (119871) : lYoon et alJ (120061) : iBrott et alJ (I2011h . 
Note, however, that fast rotation is no mandatory condition to 
achieve a nearly homogeneous stellar structure. Because the con- 
vective cores of very massive stars are very large, their remaining 
radi ative envelopes can eas ily be stripped off by strong mass loss 
(cf ., lYungelson et al. 1 120081) . E.g., for the WNh stars in our sam- 
ple we find mass-loss timescales of T w j n d = M/M ss 1..6 10 6 yr, 
similar to their evolutionary timescales of r nuc = 7.2 10 10 yr • 
(M/M O )/(L/L ) * 3.. 5 10 6 yr. In addition, other scenarios, like 
mergers due to frequent ste llar collisions in dense clusters (e.g., 
iPortegies Zwart et alJI 19991 ) might be considered in this context. 

An argument against chemical homogeneity are the low 
surface temperatures of the Arches cluster stars (30-35 kK). 
Because of their high mean molecular weight, homogeneous 
stars are expected to be very compact, and to show very high 
surface temperatures. iMartins et al.l J2009) indeed found evi- 
dence for homogeneous evolution for an extremely early WN3h 
star in the SMC, based on its high temperature (T* = 65 kK). 
This star is clearly located to the left of the main sequence, but 
has still 50% hydrogen at its surface, in agreement with evo- 



lutionary tracks for fast rotating quasi homogeneous stars by 
iMevnet & MaedeJ (120051) . So why are the Arches cluster stars 
so much cooler? 

In principle, their temperatures may indicate that the Arches 
stars are evolved ma i n sequ ence stars, and not homogeneous. 
However, llshii et all d 19991) found that very massive homo- 
geneous stars at high metallicities can have inflated stellar 
envelopes. This effect is expected to occur for stars close 
to the Eddington li mit, and to be enhanced by fast rotation 
dMaeder et al.ll2008l) . The inflation is caused by the Fe opacity 
peak at temperatures around 160 kK. The effect is thus metal- 
licity dependent. Note that the same opacity p eak is expected to 
drive the inner winds of early W R subtypes dNugis & Lamersl 
12001 IGrafener & Hamannll200l . However, for late WR sub- 
types, like the Arches stars, the opacity peak is located below 
the wind driving zone, and may thus affect the stellar envelope. 
Th e existence of the en velope extension has been questioned 
by iPetrovic et al.l (120061) because in their hydrodynamic stellar 
structure models the effect was suppressed by mass loss. It is 
thus not clear whether the envelope inflation occurs in nature, 
but it is in line with the observed low te mperatures of many H- 
free WR stars (e.g., lHamann et al.ll200"6l) . and the tendency that 
late WR subtype s are preferential ly found in high metallicity en- 
vironments (e.g.. ICrowmerj|20"07l) . 

7. Conclusions 

Based on a semi-empirical study of the most massive stars in 
the Arches cluster we could confirm key properties of WR-type 
stellar w inds that have been predicted in our previous t heoretical 
studies (IGrafener & Hamannl 12001 IVink et al.l 1201 lb . We find 
that mass loss is enhanced close to the Eddington limit, with 
a dependence of the form M(T e , T+,L, X^,Z) that shows 1) a 
strong dependence on the Eddington factor r e , 2) a weak depen- 
dence on the luminosity L, and the hydrogen surface abundance 
X£., and 3) a strong increase of mass loss for decreasing stellar 
temperatures r+. Due to the properties of the sample, the latter 
is however not fully established. The qualitative agreement be- 
tween models and observations suggests that the mass-loss prop- 
erties of WR stars are determined in analogy with the models, by 
the influence of the Fe-opacity peaks in deep atmospheric layers. 

The strong F-dependence marks an important paradigm shift 
with respect to previous, luminosity-d ependent relations, that 
are based on purely empirical studies (iNugis & Lamersl l2000fc 
lHamann et al.l 120061) . The latter are widely used, e.g., in stel- 
lar evolution models, where the WR-phase is usually identified 
by a hydrogen-deficient surface composition. According to our 
work, the physically relevant parameter, which causes WR-type 
mass loss, is the proximity to the Eddington limit. A hydrogen 
deficiency then occurs naturally as a consequence of the strong 
mass loss, potentially in combination with rotational mixing. 

As a by-product of our work we obtain information on the 
evolutionary status of the sample objects. We find strong ev- 
idence that the luminous WNh stars in the Arches cluster are 
not evolved objects, but mostly very massive, hydrogen-burning 
stars, supporting the picture that the m ost massive stars are 
WNh stars in young stellar clusters (e.g.. Ide Koter et al]|1997l : 
ICrowtheret alfcoiO) . 

The mass loss in the WR phase is of paramount impor- 
tance for the evolution, and the death of massive stars, as well 
as the formation of long GRBs. With our present work we 
have made a first step towards a more physically motivated 
description of WR-type mass loss, incorporating the depen- 
dence on the Eddington factor T e . Other predicted dependen- 
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Fig.A.l. Homogeneous stellar structure models for the mass 
range 0.6-4000 M Q (black symbols). The models are com- 
puted for hydrogen mass fractions Xu = 0.0-0.7, in steps 
of 0.1 (from top to bottom). Dashed red lines indicate the 
M-L relations for Xu = 0.1, 0.4, and 0.7 according to 
Eq.l9l&Table lA.ll and the solid red line corresponds to pure He 
models (Eq.lTUl&TablelTOV 
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Fig.A.2. Eddington parameters for the same models as in 
Fig. IA.ll The quantity G4 = T e /(1 -T e ) 4 is theoretically expected 
to follow a relation with G4 oc M 2 (see Eq. IB.5b . The symbols in 
this plot are the same as in Fig. IA.ll 

cies, e.g. on the stellar temperature, and the metallicity, demand 
for similar studies of large stellar samples in different environ- 
ments. Combined observational and theoreti cal programs, like 
the ongoing VLT-FLAM ES Tarantula Survey dEvans et al.l201 It 
iBestenlehner et al1l201 ll Bestenlehner et al. in prep.), will give 
improved insight into exactly these questions in the near future. 

Appendix A: Mass-luminosity relations for a broad 
range of stellar masses 

In Sect.|4] of this work we provided fitting formulae for 
mass-luminosity relations L(M, Xu), and the corresponding 
masses M^omiL, Xu), and Mueb(L) for massive, chemically- 
homogeneous stars. We focused on a mass range of 12-25OM , 
which covers the observed properties of the most massive stars 
known. Within this mass range we obtained a very good fit qual- 
ity for L(M, Xu), with a fitting error of up to 0.02 dex. The re- 



sulting fitting relations are however strictly limited to this mass 
range. The reason is that just within this range the M-L rela- 
tion changes from a power law L oc M n , with r\ as 3.8 for low 
masses, to a nearly linear relation with 77 — > 1 for very massive 
stars. Our fits in Sect.|4]describe this 'bending' of the M-L rela- 
tion very well, but are not capable of reproducing the asymptotic 
behaviour for higher and lower masses. We therefore decided to 
provide separate fitting relations in this appendix, that cover the 
range of 2-4000 M for H-burning stars, and 0.6-1000 M for 
He-burning stars. 

In analogy to Sect. [4] we obtain separate fits to our compu- 
tational results for H-burning stars with Xu = 0.1-0.7, and He- 
burning stars (Xu = 0). The resulting relations have the same 
form as Eqs.|9]&[T0]in Sect.|4] Moreover, we obtain separate re- 
lations for the steep low-mass part, and the flat high-mass part 
of the M-L relation. The corresponding coefficients are given in 
Table lA. ll The maximum fitting errors in log(L/L ) are 0.05 for 
H-burning stars with masses of 2-100 M Q , 0.02 for H-burning 
stars with 60-4000 M , 0.03 for He-burning stars with 0.6- 
100 M , and 0.003 for He-burning stars with 60-1000 M . The 
results are displayed in Fig. IA.ll 

Inverting these relations, we obtain Mh om (L, Xu), and 
Mueb(L) in the same way as in Sect.|4] i.e., Mh om (L,Xu) is given 
by Eqs.fTTI&fTZl and Mueb(L) by Eq.[T3] with coefficients from 
TableUO 

Appendix B: Asymptotic behaviour of the 
Eddington factor for large stellar masses 

For the highest stellar masses, the M-L relation has to become 
nearly linear, because T e — > 1 . As T e itself hardly changes in 
this range, the value of 1 — T e becomes physically more rele- 
vant. Our previous fitting relations, however, chiefly focus on 
M, and L, i.e., on the correct reproduction of F e . The resulting 
value of 1 - T e can be significantly affected by small uncertain- 
ties due to fitting errors. E.g., for the largest masses of up to 
4000 M , our relation in Eq.|9]reproduces the stellar luminosity 
and thus F e rather precisely, but the error in 1 — T e amounts to 
0. 1 dex. To improve on this we take advantage of a theoretically 
predicted scaling relation th at has recently been pointed out by 
lOwocki & Shaviv| (|in press!) , and goes back to the original work 
bv lEddingtonl (119181) . 

With Eddington's original definition 

Q-& S? ? = ^F> ( B -i) 

the total pressure P due to gas + radiation can be written in the 
form 

Assuming/? is a constant, this relation leads to a poly tropic equa- 
tion of state with P oc p 4 / 3 , which is the basis of the Eddington 
Standard Model. 

Eddington's assumption that f3 is a constant, holds very well 
in large parts of the stellar envelope where the energy transport is 
dominated by radiation, and where the opacity is dominated by 
the constant electron scattering opacity Xe- In that case it follows 
from the equation of radiative diffusion that 

dP ra d XzpL r GM(r)p dP 
or 4nr~c r L or 
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Table A.l. Mass-luminosity relations for chemically-homogeneous stars. 



No. 


ev. ph. 


Eq. 


mass range 


Fi 


F 2 


F 3 


^4 


F s 


F 6 


Fi 


^8 


F 9 


1 


H-b. 


121 


12-250 Mq 


2.875 


-3.966 


2.496 


2.652 


-0.310 


-0.511 








2 


H-b. 


® 


2-100 M 


1.967 


-2.943 


3.755 


1.206 


-0.727 


-0.026 








3 


H-b. 


El 


60-4000 Mq 


3.862 


-2.486 


1.527 


1.247 


-0.076 


-0.183 








4 


H-b. 


IB~7l 


2-30 M Q 


-2.688 


-7.843 


2.471 


2.758 


-0.233 


-0.747 








5 


H-b. 


lB~7l 


12^1000 Mq 


-2.416 


-5.118 


1.869 


-0.400 


0.064 


0.050 








6 


He-b. 


[TOl 


8-250 M Q 


3.017 


2.446 


-0.306 














7 


He-b. 


[Tol 


0.6-100 M B 


2.635 


2.986 


-0.488 














8 


He-b. 


[TOl 


60-1000 M Q 


3.826 


1.619 


-0.099 














9 


He-b. 


lB~8l 


0.3-100 M e 


-2.204 


1.831 


0.149 














10 


He-b. 


IB~8l 


12-500 Mq 


-1.676 


1.075 


0.404 














11 


H-b. 


edeh 


12-250 Mq 


4.026 


4.277 


-1.0 


25.48 


36.93 


-2.792 


-3.226 


-5.317 


1.648 


12 


H-b. 


E0E3 


2-100 Mq 


2.582 


0.829 


-1.0 


9.375 


0.333 


0.543 


-1.376 


-0.049 


0.036 


13 


H-b. 


Eiiin] 


60-4000 Mq 


10.05 


8.204 


-1.0 


151.7 


254.5 


-11.46 


-13.16 


-31.68 


2.408 


14 


H-b. 


IB.13IIB.14I 


2-30 Mq 


5.303 


5.918 


-1.0 


16.58 


-4.292 


-72.89 


-7.881 


-13.76 


3.206 


15 


H-b. 


IB.13IIB.14I 


12^1000 Mq 


-14.60 


3.125 


1.0 


251.0 


15.63 


72.24 


18.20 


12.21 


0.781 


16 


He-b. 


EH 


8—250 Mq 


3.997 


-1.0 


25.83 


-3.268 












17 


He-b. 


ED 


0.6—100 Mq 


3.059 


-1.0 


14.76 


-2.049 












18 


He-b. 


EH 


60-1000 Mq 


8.177 


-1.0 


105.5 


-10.10 












19 


He-b. 


|B.15| 


0.3-100 Mq 


-6.144 


1.0 


52.54 


6.711 












20 


He-b. 


IB.15I 


12-500 Mq 


-1.330 


1.0 


5.919 


2.475 













Notes. The table describes the mass-luminosity relations L(M, X a ), and M(L, X a ) that are derived in this work, as well as the relations T C (M, X H ), 
and M(r e ,Xn) in the appendix. For each relation we give the assumed evolutionary phase, the corresponding equation numbers, the relevant mass 
range, and the numeric coefficients F„ . 



The latter equality follows from the definition of the classical 
Eddington factor F e (Eq.Q}, and the equation of hydrostatic equi- 
librium. Under the assumption that T e is a constant, Eq lB.3l is just 
the derivative of the relation P la d = T e P, i.e., it becomes equiv- 
alent to Eq.|BT] We thus have P md /P = (1 - j3) = F e , and with 
EqE2] 



'39* 



(B.4) 



Here we are interested in how F e changes qualitatively with 
the stellar mass M. For this purpose we express P, and p in 
Eq. IB.4l bv two-point scaling relations in M, and R. From hydro- 
static equilibrium it follows directly that P/R * GMp/R 2 . With 
p =a M/R 3 this leads to P « M 2 /R 4 . Using the last two relations 
with Eq. IB.41 we obtain the approximate scaling relation 



(B.5) 



1, r r. 

oc — oc M . 

d-r e ) 4 p 4 



Based on this simple scaling law we obtain very precise fit- 
ting relations for the quantity 



G 4 = 



d-r e ) 4 



(B.6) 



The resulting fits to our computational results are plotted in 
Fig. IA.2l For H-burning, chemically-homogeneous stars with 
Xu = 0.1-0.7 we use a relation of the form 



log(G 4 ) = [Fi + F 2 log(l 
+ [F 3 + F 4 log(l 
+ [F 5 +F 6 log(l 



■X H )] log(M/M Q ) 
■X H )] log(M/M ) 2 



(B.7) 



In Table lA. ll we give the coefficients F\-F(„ for the mass ranges 
2-30 M Q , and 1 2-4000 M . For pure He stars (X H = 0) we use 



log(G 4 ) = Fi+F 2 log(M) + F 3 log(M) 2 , 



(B.8) 



with coefficients F\-F 3 from Table lA. 11 for the mass ranges 
of 0.3-100 M Q , and 12-500 M & . The maximum fitting error in 
log(G 4 ) for relations $BJ]&\BM is 0.03. 

For given G4, F e can be obtained from the solution of Eq.| 



1 



C 
15m 



3 VdgT 



1 - 1 



with the definitions 

C = VI7I2, 



D = (C/3 V256G 4 + 27 + 1/2) /G, 
E = 3G 4 D 2/3 -4. 



(B.9) 



(B.10) 
(B.ll) 
(B.12) 



Finally, we obtain relations for the corresponding masses 
M hom (G 4 ,X H ), and M He b(G 4 ) from Eqs.E3&EH 



log(M hom /M ) = 



with 



Fi +F 2 \og(l+X H ) + F 3 ^ 
1 +F 9 log(l +X H ) 



g = F 4 + F 5 log(G 4 ) + F 6 log(l + X H ) 2 
+ CF 7 + J F 8 log(G 4 ))log(l +X H ) 

and 



(B.13) 



(B.14) 



(B.15) 



log(MHeb/M ) = F 1+ F 2 ^F i + F 4 \og(G 4 ). 
The coefficients F„ are given in Table lA. ll 
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